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We present a quantitative analysis of the microscopic Dirac spectrum which is complex in the presence of a non- 
vanishing quark chemical potential. Data from quenched SU(3) lattice simulations for different volumes V and 
small values of the chemical potential ^ are compared to analytical predictions from matrix models. We confirm 
the existence of two distinct limits for weakly and strongly nonhermitian Dirac operators. Good agreement is 
found in both limits, confirming the different scaling of chemical potential and eigenvalues with the volume. 



1. Introduction 

Lattice QCD at nonzero density remains a diffi- 
cult topic. Several methods to solve the problems 
related to the complex phase of the Dirac oper- 
ator determinant have recently been suggested: 
multiparameter reweighting, Taylor expansion, 
and analytic continuation from imaginary /i 
So far, all three methods work only in the vicinity 
of the phase transition close to the temperature- 
axis. Here, we instead explore the /i-axis at zero 
temperature and small fi. In this region matrix 
models analytically predict the complex micro- 
scopic Dirac spectrum [2]. We test these predic- 
tions for pure SU(3) gauge theory [3], where ^ ^ 
is easily implemented. Although the quenched 
approximation is problematic 0], our test is an 
important check, as the matrix model remains 
valid and predictive for Nf > flavors. For exam- 
ple, the 1-flavor matrix model has been used to 
test a new factorization algorithm for the com- 
plex action problem |^. In fact, at fi — the 
situation is analogous: while most tests of matrix 
model predictions have been performed only for 
pure gauge theories, dynamical fermions are also 
well described [H] ■ For a review of matrix models 
in QCD we refer to Ref. [7]. 



2. Matrix model predictions 

In Ref. PI a matrix model has been formulated 
in terms of N complex eigenvalues and solved in 
the limit iV oc U — > oo. It has the same global 
symmetries as QCD at /i 7^ due to an equiva- 
lence to the matrix model of Ref. ^4; at small 
^ in the phase with broken chiral symmetry. We 
give here only the quenched results for the mi- 
croscopic density and refer to Ref. El for more 
details and higher-order correlation functions. 

In Ref. [2] two different large-TV limits were 
found, those of weak and of strong nonhermiticity. 
In QCD at ^ = 0, the lowest Dirac eigenvalues 
scale with 1/V to build up a finite condensate 
according to the Banks-Casher relation. In the 
weak- nonhermiticity limit, first introduced in [S], 
this scaling remains unchanged. The support of 
the density remains quasi one-dimensional as we 
send ^ ^ 0, keeping limAr^oo 2A^/x^ = fixed. 
In the strong- nonhermiticity limit, the eigenval- 
ues fill a two-dimensional surface, and thus the 
scaling is modified to l/VV. We find for the 
quenched microscopic density at weak nonher- 
miticity in the sector of topological charge i/ 
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Here, we have rescaled the complex eigenvalues 
according to ^ ~ Vz. At strong nonhermiticity, 
fi is kept fixed as V oo, and one obtains [2] 



Pst 



rong 



(0 = 



ICI exp 



2^2 



1^ 

2/i2 



.(2) 



rescaling the eigenvalues according to ^ ^ yVz. 
The microscopic density Eq. is rotationally 
invariant in contrast to the macroscopic support. 

3. Lattice Data at /i 7^ 

Our lattice calculations were done using the 
staggered Dirac operator with chemical potential 
/i 7^ for a pure SU(3) gauge theory, correspond- 
ing to Nf — 0. There are two reasons for this 
choice. First, we need high statistics (w 20, 000 
configurations^) to test the matrix model predic- 
tions of Eqs. 1)112(1 . Ginsparg- Wilson type opera- 
tors would be too expensive here. Second, Wil- 
son fermions break chiral symmetry explicitly and 
have complex eigenvalues already at fi ~ 0. 

In the simulations, we have chosen (3 = Qj ~ 
5.0 in the strong-coupling regime for the follow- 
ing reason. At zero temperature and /i = the 
matrix model is equivalent to a low-energy 
effective theory of QCD (the effective chiral La- 
grangian) in the so-called e-regime, where the 
zero-momentum modes of the pseudo-Goldstone 
bosons dominate (see e.g. 171). When the nonzero- 
momentum modes start contributing, this equiv- 
alence breaks down. The corresponding scale, the 
so-called Thouless energy, is a function of both (3 
and V . As /? is increased, fewer eigenvalues are 
described by the matrix model, an effect that can 
be compensated by increasing V . For our small 
volumes V = 6^, 8^, 10'*, however, we are limited 
to relatively small values of (3 (assuming a similar 
behavior of the Thouless energy for small /i 7^ 0) . 

At strong coupling and thus away from the con- 
tinuum limit, staggered fermions have the disad- 
vantage of shifting the (topological) Dirac zero 
modes and mixing them with the nonzero modes. 
We have accounted for this by setting v = {] in 
Eqs. (QJ and 10) above. 
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iThe y = 10** data in Figs. EEl contain only 4,000 config 
urations, therefore the statistical errors are larger. 



Figure 1. Densities of small Dirac eigenvalues cut 
along the real axis (left) and parallel to the imag- 
inary axis at the first maximum (right). Data 
(histograms) for y = 6'* at ^ = 0.006 (top), 
T/ = 8* at ^ = 0.003375 (middle), and V = 10"^ 
at /i = 0.00216 (bottom) vs. Eq. (QJ for weak 
nonhermiticity. 



The three different lattice volumes are cho- 
sen in order to test the different scaling behavior 
of eigenvalues and chemical potential for weakly 
and strongly nonhermitian lattice Dirac opera- 
tors. For more simulation details we refer to 0- 

Let us first discuss the data at weak nonher- 
miticity. In Eq. the average level spacing be- 
tween consecutive eigenvalues is tt. In order to 
compare with lattice data, we first determine the 
mean level spacing d from the data averaged over 
many configurations and then rescale the lattice 
eigenvalues z by ^ = -nzj d. At the same time the 
spacing d oc 1/V provides us with the weak non- 
hermiticity parameter a = fi y^2/d ~ 0.19 to be 
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Figure 2. Densities of small Dirac eigenvalues cut 
along the real (left) and imaginary (right) axes. 
Data (histograms) for V — 6* (top), V — 8* (mid- 
dle), and V = 10"* (bottom), all at ^ = 0.2, vs. 
Eq. Q for strong nonhermiticity. 



used in Eq. After normalizing the histograms 
to unity we thus obtain the parameter-free com- 
parison shown in Fig. The three different val- 
ues of n are such that the product Vn^ is con- 
stant. All data are well described by the same 
value^ of a determined from the volume V = 6'^. 
Apart from the scaling Vfi^ oc o? we have thus 
also confirmed the scaling of the complex Dirac 
eigenvalues with the volume. 

At strong nonhermiticity we have chosen [i = 
0.2 for all lattice volumes. In this case, the eigen- 
values are rescaled with the square root of the vol- 
ume, ^ = TTz/^fd, where the level spacing is now 



^The values of a obtained from the V ■ 
agree within errors. 



8* and lO"* data 



obtained by averaging over the smallest geomet- 
ric distance between complex eigenvalues. The 
normalized data confirm the prediction of Eq. |(2Jl 
together with the different scaling. 

In conclusion, we have confirmed the analytical 
matrix model predictions in two different scaling 
regimes from quenched QCD lattice data. For 
the small values of /i at weak nonhermiticity, un- 
quenched simulations should be feasible as well, 
keeping V jj? fixed. On the other hand, the lat- 
tice simulations m close to the phase transition 
are at strong nonhermiticity, since /i is kept con- 
stant. So far, our analytical predictions are re- 
stricted to the broken phase and, unfortunately, 
are not suitable to describe the transition region. 
It would be very interesting to compute predic- 
tions for the unbroken phase as well. 
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